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ABSTRACT 

We report the discovery of the first hyperluminous infrared galaxy (HyLIG) in the course 
of the European Large Area ISO Survey (ELAIS). This object has been detected by ISO 
at 6.7, 15, and 90 ^im, and is found to be a broad-line, radio-quiet quasar at a redshift: z 
= 1.099. From a detailed multi-component model fit of the spectral energy distribution, we 
derive a total IR luminosity: Lir (1-1000 ^m) « 1.0 x 10^^ hg^ Lq (qg = 0.5), and discuss 
the possible existence of a starburst contributing to the far-IR output. Observations to date 
present no evidence for lens magnification. This galaxy is one of the very few HyLIGs with 
an X-ray detection. On the basis of its soft X-ray properties, we suggest that this broad-line 
object may be the face-on analogue of narrow-line, Seyfert-like HyLIGs. 

Key words: Galaxies: Infrared — Galaxies: Quasars: Spectral Energy Distri- 
bution — Galaxies: Surveys — Galaxies: Individual ([CCS88] 163831.4-f411107; 
ELAISP90_J164010-f410502) — Dust 
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Sanders & Mirabel 1996 for a review and definition of the 
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Universe (Sanders & Mirabel 1996 and references therein) 
but also at moderately high redshift (Small et al. 1998). 
While a contribution to the IR output can arise both from 
a dust-enshrouded active galactic nucleus (AGN) and from 
a (generally circumnuclear) starburst, there is strong evi- 
dence from optical or mid-IR spectroscopy for an increasing 
fraction of objects with AGN-like characteristics when pro- 
gressing toward higher luminosities (e.g., Lutz et al. 1998; 
Veilleux, Sanders, & Kim 1999). 

Of particular relevance for a better understanding of 
this important population are the objects at the upper end 
of the luminosity distribution (the so-called hyperluminous 

infrared galaxies; hereafter HyLIGs).P] If their luminosities 
are starburst-dominated, they are potentially forming the 
bulk of their stellar population in a single, violent episode. 
This makes them an interesting sub-population of galax- 
ies, but work so far is severely hampered by inhomogenous 
and/or AGN-biased selection. Only 13 HyLIGs are known 
at present from far-IR or sub-mm surveys, and a further 
12 from cross-correlating far-IR source lists with AGN cat- 
alogues (Rowan- Robinson 2000; hereafter RR). 

In regard of the exceedingly small space density of 
HyLIGs (van der Werf et al. 1999; RR), only major ex- 
tragalactic IR surveys offer good prospects to uncover new 
candidates. Of particular interest in this respect is the Euro- 
pean Large Area 750 Survey (ELAIS) which was conducted 
over ~ 11 square degrees at 15 fim and 90 ^m (down to 
~ 3 mjy and ~ 100 mjy, respectively) and over 6 square 
degrees at 6.7 ^m (down to « 1 mJy). We refer the reader 
to Oliver et al. (2000; Paper I) for a complete description of 
this project. Details on the technical aspects of the CAM and 
PHOT observations can be found in Serjeant et al. (2000; 
Paper II) and Efstathiou et al. (2000a; Paper HI), respec- 
tively. Assuming pure luminosity evolution of the form: (1 + 
zf-^ (Boyle, Shanks, & Peterson 1988; Saunders et al. 1990), 
we estimate that about 3 HyLIGs might be detected in the 
total area covered by ELAIS. This estimate is considerably 
uncertain, as it is strongly dependent on the spectral energy 
distribution (SED) and evolutionary model assumed. As we 
can safely rule out no-evolution models, however, we find 
that the formal significance of detecting one object is very 
high (« 35a). We report in this paper on the discovery of 
such a candidate HyLIG in the ELAIS N2 region. 



2 IDENTIFICATION OF THE NEW HYLIG 

This new HyLIG was discovered as part of a study aiming at 
studying the IR colour properties of ELAIS galaxies (Morel 
et al. 2001). The procedure used to isolate galaxies poten- 
tially detected at 6.7, 15, and 90 fim was to cross-correlate 
the CAM and PHOT source lists using a 35" association 
radius. Details on the production of the CAM and PHOT 
catalogues can be found in Paper II and Heraudeau et al. 
(2001), respectively. 

Figure 1 shows the PHOT error circle of the HyLIG 
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Figure 1. Optical (r') contour map of ELAIS J1640-I-41 centered 
on the PHOT detection (cross). The lowest contour is drawn at 
3cr, and then by steps of 6cr. The error circle for PHOT (35"), as 
well as the 6" error circles for the detections at 6.7 fim (solid) and 
15 fim (dashed) are also overlayed. The new HyLIG is indicated 
by an arrow. 



(ELAISP90_J164010-|-410502; hereafter ELAIS J1640+41) 
overlayed on a deep (down to r' ~ 24) optical image. The 
CAM error circles are also shown. The only obvious optical 
counterpart to the ISO sources is found to be a point-like 
object at: a = 16''40'"10.16" and S = -h41°05'22.3" (J2000). 
This source is, within the uncertainties, spatially coincident 
with an optically-selected quasar at a quoted redshift of 
1.097 ([CCS88] 163831.4-^411107; Crampton et al. 1988). 
As can be seen in Fig.l, a bright optical source (as well as a 
number of much fainter ones) appear to lie in close vicinity 
of the quasar. Since no redshifts are available for these ob- 
jects, it is unclear at this stage whether they are physically 
associated with the new HyLIG. 

We used the quasar B-band source counts of Wisotzki et 
al. (2000) to estimate the number of random PHOT-quasar 
associations in our sample at B = 17.2 mag (see Table 2), or 
brighter. The expected total number of random associations 
in the search radius is only 0.011 for the 285 PHOT sources. 
Furthermore, 35% of the ELAIS N2 area is covered by the 
parent sample of the quasar, i.e., only about 10% of the total 
ELAIS area. The total number of random associations in the 
overlap region is therefore only « 1 x 10~^, giving strong 
support for the reality of the PHOT-quasar association. 

The quasar does not appear in the IRAS FSC or FSR 
catalogues. To test the reliability of our PHOT detection, 
we made IRAS ADDSCANs of the source, as well as of 15 
galaxies with reliable 90 /im detections at similar fiux densi- 
ties (60-80 mJy) from the Paper HI catalogue, and 16 control 
positions at random 20' offsets from the target positions. We 
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Table 1. Line properties of ELAIS J1640+41 



Table 2. Photometric properties of ELAIS J1640+41 



Line 



Relative line flux 



Wrest 

(A) 



FWHM 
(km s"-"-) 



Band 



Flux density/Mag. Source/Ref. 



C III] AigOQ" 8.0 ± 4.0 19.7 ± 3.0 6800 ± 300 

C II] A2327 1.0 4.40 ± 2.0 3500 ± 600 

Mg II A2798 3.5 ± 1.5 21.0 ± 3.0 3000 ± 300 

NOTES TO TABLE 1 — 

°: Deblended from Al III A1857. Might receive a contribution 

from Si III] A1892. 




2000 aSOO 3000 3500 4000 

ReHl— fr'tirne wtiVelengLh (A) 

Figure 2. Rest-frame UV/optical spectrum of ELAIS J1640-I-41 
(assuming a redshift of 1.099). The main lines are indicated (the ® 
symbols mark terrestrial telluric absorption features). The spec- 
trum has been corrected for both atmospheric and foreground 
galactic extinction. We assumed E(B — V) = 0.006 mag (Schlegel, 
Finkbeiner, & Davis 1998). 



found our PHOT 90 /im data to be not well reproduced by 
the ADDSCAN 100 iJ,ra measurements, although most of the 
flux densities were positive. However, the control fields gave 
values distributed around zero as expected, from which we 
derived a. ^ 4a ADDSCAN detection of our PHOT source. 
This confirms the reality of the source, but we also conclude 
that it is too faint to reliably measure the flux density with 
ADDSCANs. We only consider in the following upper limits 
from the FSC. 

Optical spectroscopy was obtained at the Nordic Opti- 
cal Telescope using ALFOSCn {R ~ 700) in non-photometric 
conditions and seeing of 1.5' . We used a 1.2"-wide slit ori- 
ented EW. The exposure time was 3 x 900 s. The spectrum 
was reduced using standard IRAF|S| routines and is shown in 
Figure 2. We derive a redshift: z = 1.099 ± 0.002. This value 
is consistent with previous estimates (Crampton, Cowley, & 
Hartwick 1989), and is used in the following. Table 1 gives 
the relative line fluxes, EWs, and FWHMs. 



■f The Nordic Optical Telescope (NOT) is operated on the island 
of La Palma jointly by Denmark, Finland, Iceland, Norway, and 
Sweden, in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofi'sica de Canarias. ALFOSC is owned by 
the Instituto de Astrofi'sica de Andalucia (lAA) and operated at 
the NOT under agreement between lAA and the NBIfAFG of the 
Astronomical Observatory of Copenhagen. 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the Na- 
tional Science Foundation. 



0.1-2.4 keV 

U 

r' 
R 
J 
H 

Ks 

ISOCAM-LW2'= 
12 fim 

ISOCAM-LW3= 
25 /jm 
60 /im 

ISOPHOT-C90 
100 ^lra 
450 fira 
850 ^J.ra 
4.85 GHz 
1.4 GHz 



1.9 X 10-^5 mjy 
16.362 ± 0.200 
17.17 ± 0.25 
16.94 ± 0.05 
16.686 ± 0.200 
16.395 ± 0.134 
15.520 ± 0.139 
15.164 ± 0.180 
4.99 ± 0.66 mJy 
<200 mJy {3a) 
7.75 ± 0.72 mJy 
<200 mJy (3o-) 
<200 mJy {3a) 
72 ± 23 mJy 
<1000 mJy (3a) 
<81.6 mJy (3a) 
<6.96 mJy (3a) 
<18 mJy {5a) 
<0.95 mJy (5a) 
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Morel et al. (2001) 
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IRAS FSC 
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NOTES TO TABLE 2 — 

The optical and near-IR magnitudes have been corrected for 
foreground galactic extinction assuming E{B — V) = 0.006 mag 
(Schlegel et al. 1998) and the conversion factors of Cardelli, Clay- 
ton, & Mathis (1989). 

'^: Assuming continuum emission with a photon index F = 2 (see 
Reeves & Turner 2000) and a foreground galactic extinction with 
a column density: A^h = 1-1 X 10^° cm-^ (Stark et al. 1992). 
^■. We do not adopt the value quoted by Crampton et al. (1988): B 
= 17.9 ± 0.2 mag, as it is discrepant both with APM and USNO 
magnitudes. This might be partly due to quasar variability. 
'^: These values are subject to a possible systematic scaling by up 
to a factor of about 1.5 (see Paper II). 



3 MODEL FIT OF THE SED 

Table 2 summarizes the available photometric data for this 
seldom-studied object. New sub-mm observations (see Far- 
rah et al. 2001) were carried out with the Submillimetre 
Common-User Bolometer Array (SCUBA) at the James 
Clerk Maxwell Telescope (JCMT) in photometry mode. 
Data reduction was performed using the standard SURF 
pipeline software. The source is undetected both at 450 and 
850 pim, with 3cr upper limits of 81.6 and 6.96 mJy, respec- 
tively. The region surrounding the HyLIG has unfortunately 
not been observed as part of the radio survey of the north- 
ern ELAIS areas (Ciliegi et al. 1999), but the available ra- 
dio flux limits imply that this quasar is radio-quiet (Sopp 
& Alexander 1991). Following Veron-Cetty & Veron (2000), 
we obtain: Mb ~ - 26.0 mag. 

The SED is displayed in Figure 3. Assuming that the 
host galaxy is similar to coeval quasars (e.g., McLure et 
al. 1999), we used the radio galaxy K-z relation of Eales 
et al. (1997) to estimate the host galaxy near-IR magni- 
tude, obtaining: K — 17.0 ± 0.4 mag. As can be seen in 
Fig. 3, the contribution of the host galaxy to the SED is 
negligible. In order to estimate the energy budgets to the 
IR luminosity, we modelled the SED by simultaneously con- 
sidering models for dust-enshrouded AGNs (Efstathiou & 
Rowan-Robinson 1995) and starbursts (Efstathiou, Rowan- 
Robinson, & Siebenmorgen 2000b). Although we do not 
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Figure 3. Rest-frame SED of ELAIS J1640-I-41 from X-rays to 
radio (data from Table 2). Tlie SED is fitted with a combination 
of one of the starburst models of Efstathiou et al. (2000b) {dashed 
line; ry = 50, t = 26 Myrs) and an AGN dusty torus model [dotted 
line) from Efstathiou & Rowan-Robinson (1995). We assume a 
half-opening angle for the torus of 45° and a face-on inclination. 
The starburst SED is extrapolated in the radio using the far- 
IR/radio correlation (Helou, Soifer, & Rowan-Robinson 1985). 
For the host galaxy, the SED of a 4 Gyr old stellar population 
{dotted- dashed line; Fioc & Rocca-Volmerange 1997) has been 
normalized to K = 17.0 ± 0.4 mag (see text). 



Figure 4. Comparison of the rest-frame SEDs of a number of 
HyLIGs with optical evidence for quasar activity (data from RR 
and references therein). The solid and dashed lines show the 
model fitted to ELAIS J1640-I-41 (see Fig.3) and the mean SED 
of radio-quiet quasars in the UVSX sample (Elvis et al. 1994), 
respectively. The data points are normalized to their rest-frame 
flux at 15 fim. The vertical scale is arbitrary. 



Mg II absorbers in high-quality, moderate-resolution spec- 
troscopic data (e.g., Goodrich et al. 1996). 



claim to have exhaustively explored the possible physical 
parameter space of these models, we found that the SED 
is most naturally explained by a combination of these two 
components. Pure AGN models have difficulties in account- 
ing for the far-IR emission. We achieve a satisfactory fit to 
the SED (see Fig.3) by assuming roughly similar contribu- 
tions from dust grains heated by the nuclear UV/optical 
continuum (Ljj^^'^ ~ 5.2 x lO^'^ hgc, Lq) and by a starburst 
{L^i ^ 5.0 X 10^^ h^^ Lq). The difliculties in fitting the 
region around 3-7 ^ra might be related to uncertainties in 
the ISOCAM absolute calibration (Paper II). 

The total rest-frame IR luminosity we derive: Lir (1- 
1000 lira) « 1.02 x lO" hg^ L©, is slightly above the 
threshold for a HyLIG classification in the scheme of RR. 
The robustness of the derived luminosity has been explored 
by considering two low-redshift objects in the list of RR 
with detections in all IRAS bands (IRAS 07380-2342 and 
IRAS 18216-1-6418). Similarly to the approach adopted here, 
the IR luminosities have been derived by RR from a multi- 
component model fit of the SED, and are found to agree 
to within 25% with the values determined by the method 
prescribed by Sanders & Mirabel (1996). We conclude that 
estimates of the total IR luminosity are fairly insensitive to 
the choice of the method used, although the relatively low 
value we obtain for ELAIS J1640-I-41 (coupled with the quite 
large uncertainties in the IR fluxes) indicates that this ob- 
ject can only be formally considered as a candidate HyLIG 
at this stage. 

Gravitational lensing is known to enhance the luminos- 
ity in a number of HyLIGs (e.g., Serjeant et al. 1995). While 
our data do not allow us to rule out lens magniflcation at 
this stage, this possibility can be explored in the future via 
high-resolution HST imaging or by seeking for intervening 



4 DISCUSSION 

4.1 ELAIS J1640+41 in relation to other HyLIGs 

We compare in Figure 4 the SED of ELAIS J1640-I-41 with 
a sample of HyLIGs with optical evidence for quasar ac- 
tivity. It can be seen that ELAIS J1640-I-41 (which is the 
least luminous in this subset) presents a relative excess in 
the far-IR. Interestingly, the strongest outlier (BR 1202- 
0725) is also the most luminous (RR) . This trend of hotter 
colour temperature with increasing luminosity might sug- 
gest an increasing contribution from AGN emission at high 
luminosities (see also Haas et al. 2000). However, such an 
effect is not apparent in the sample of quasars studied by 
PoUetta et al. (2000). The detection of a large reservoir of 
molecular gas in BR 1202-0725 (Ohta et al. 1996), but not 
in PG 1634-1-706 (Barvainis et al. 1998) also seems to be at 
variance with this simple picture. 

We also show in Fig. 4 the mean SED for radio-quiet 
quasars in the UVSX sample (Elvis et al. 1994). It can be 
seen that ELAIS J1640-I-41 displays a significant far-IR ex- 
cess compared to this population, hence giving credence to 
the existence of a starburst contributing to the IR output. 
Indicative of the unsually strong level of IR emission from 
ELAIS J 1640-1-41 is the fact that a number of more pow- 
erful quasars in the optical regime are not member of the 
HyLIG class (McMahon et al. 1999). Assuming a Salpeter 
initial mass function, we derive from the luminosity of the 
starburst component a current star-formation rate in the 
range: 1-7 x 10^ h^^ M© yr~^, depending upon the stellar 
low- and high-mass cutoffs assumed (see Thronson & Telesco 
1986). Although such high values have already been claimed 
in HyLIGs with quasar activity (e.g., McMahon et al. 1999), 
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we caution that the paucity of far-IR detections does not al- 
low us to set stringent constraints on the star-formation rate 
in ELAIS J1640-H41. 

Current samples of HyLIGs are severely affected by se- 
lection effects. General conclusions regarding the character- 
istics of this population must thus await the completion of 
large, unbiased surveys. With this limitation in mind, there 
is some indication at present for a diversity of power sup- 
plies. While AGN-related emission is able to account for 
the IR output in some HyLIGs (e.g., Granato, Danese, & 
Franceschini 1996), others show evidence for a dominant 
contribution from star formation (e.g., Ivison et al. 1998). 
The results of our two-component model fit to the SED of 
ELAIS J1640+41 supports a composite nature for this ob- 
ject. The existence of a number of HyLIGs whose IR lumi- 
nosity can be accounted for by a combination of an AGN 
and a starburst component has been postulated on similar 
grounds by RR and Verma et al. (2001). The inference of a 
starburst component mostly contributing to the rest-frame 
flux longward of about 50 /im is independently supported 
in several cases by the large amount of molecular material 
deduced from radio CO observations (RR). 

4.2 Soft X-ray properties of ELAIS J1640-|-41 

Inspection of the ROSAT all-sky survey 0.1-2.4 keV photon 
maps shows that our object lies within the 90% error box 
of an X-ray source detected at the 3.4a" confidence level. 
We obtain a rest-frame (uncorrected for intrinsic absorp- 
tion) X-ray luminosity in the 0.1-2.4 keV band: Lx ~ 3.7 
X 10** ergs s~^, which translates into a Lx/iboi ratio of 
about 9.4 X 10"'^. This ratio supports an AGN interpreta- 
tion for the soft X-ray emission (Wilman et al. 1998), but 
is almost two orders of magnitude higher than the upper 
limits found for HyLIGs with Seyfert 2-like optical spectra 
(Lawrence et al. 1994; Fabian et al. 1996; Wilman et al. 
1998). The spectral index of the UV/optical power-law con- 
tinuum in ELAIS J1640-I-41 {a ^ - 0.04; /^ oc i/") is typical 
of UV-selected quasars (e.g., Natali et al. 1998), and sug- 
gests little dust obscuration along our line-of-sight to the 
nucleus. In contrast, we estimate that column densities: Mn 
^ 2 X 10^^ cm~^ would attenuate the (redshifted) « 3-keV 
emission of ELAIS J1640-I-41 to levels comparable to what 
is observed in Seyfert 2-like HyLIGs.pi This value compares 
well with the typical amount of obscuring material intrinsic 
to these objects (e.g., Wilman et al. 1998). In the context of 
the unification scheme, the soft X-ray properties of ELAIS 
J1640-I-41 are thus consistent with this galaxy being a face- 
on analogue of such narrow-line HyLIGs (see Franceschini 
et al. 2000). This interpretation is independently supported 
in several cases by the detection of strongly polarized con- 
tinuum emission in edge-on objects (e.g., Hines et al. 1995). 
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